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Abstract 

A simplified Walecka-type model is investigated in a cosmological scenario. The model includes fermionic, 
scalar and vector fields as sources. It is shown that their interactions, taking place in a Robertson- Walker 
metric, could be responsible for the transition of accelerated-decelerated periods in the early universe and a 
current accelerated regime. It is also discussed the role of the fermionic field as the promoter of the accelerated 
regimes in the early and the late stages of the universe. 

1 Introduction 

One of the most fundamental discoveries in observational cosmology showed that the universe is not only expanding 
but it does at an accelerated rate [Tj. As the source of this expanded acceleration remains unknown, several 
candidates have been proposed to account for it. They include, among others, cosmological constant, matter with 
exotic equation of state, scalar fields [2], as an attempt to describe the early stages of the inflationary universe and 
current acceleration. Alternatively, fermionic fields [3j [4] have been successfully used, leading to earlier accelerated 
era as well as later decelerated and accelerated stages. These fermionic sources have been investigated by using 
several approaches, including numerical and exact solutions, perturbations, dark spinors, anisotropy-to-isotropy 
scenarios and cyclic cosmologies (see, for example [31 H])- When considering these models, a key point is the choice 
of the interaction potentials. In previous works several self-interaction fermionic potentials were tested [3J 21 [SJ |5] , 
like Nambu-Jona-Lasinio, and Yukawa type [7j. In [6] the authors developed a model in which a fermionic field 
interacting through an Yukawa-type potential led to different dynamical regimes of the universe. Another possibility 
that replaces dark energy makes use of noninteracting massless vector fields alone to derive inflation [5]. Our 
approach to cosmological dynamics is inspired by a relativistic mean field model known as Walecka model 
In the Walecka model, the relativistic nucleus interactions occur via the exchange of virtual mesons, with the 
Lagrangian controlling these interactions with additional terms [9]. The nucleons obey the Dirac equation, while 
the (scalar) virtual mesons obey the Klein-Gordon equation. Although the theory structure is invoked here, we 
suggest a completely different interpretation, reinforced by the fact that we are supposing an early universe, leaving 
the inflationary period. The consequences of these interactions to the fate of the universe evolution is one point of 
discussion in this work. Hence, our goal is to investigate whether the presence of a massless fermionic, a massive 
scalar and a massive vector fields interacting in the Robertson- Walker metric may have led to different cosmological 
regimes. Our model shows that this is indeed the case in that all three regimes - initial acceleration, deceleration 
and later acceleration - are present. One important consequence is that the fermionic field is the promoter of the 
accelerated regimes in the early and the late stages of the universe. 



2 The Walecka-type model 

We consider that the sources of the gravitational field are related to: 
(a) the Lagrangian density of a massless fermionic field with self interaction potential V(ipip) 

C f = 1 -®T»D^ - (D^)T^] - V(W>), (1) 
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where -0 and ip = ip^j represent the spinor field and its adjoint, respectively. Due to the principle of general 
covariance, the Dirac-Pauli matrices 7° are replaced by their generalized versions = e^ a , where are 
the tetrad fields, and the T M matrices satisfy the generalized Clifford algebra T"} = 2g fil/ '. Following the 
principle of general covariance, the Latin index corresponds to the local Lorentz frame and the Greek one to 
the general frame. Furthermore, 7° are the 4x4 Pauli-Dirac matrices 

where a 1 are the 2x2 Pauli matrices. The covariant derivatives in ([!} read 

= d^ip - + iqAf.tp, (3) 
Df/tp = d,4> + ^fjClfi - tq^4„ . (4) 



Above q is a constant which couples the fermionic field with the vector field A^. Moreover, f2 M is the spin 
connection 

n, = -\g pa [T^-e p b {d,e b & )]T s T\ (5) 
with T^ x denoting the Christoffcl symbols; 

(b) the Lagrangian density of a massive scalar field (f> without self interaction potential 

C b = ^d^d^-\nl4> 2 , (6) 

with m\y denoting the mass of the scalar field; 

(c) the Lagrangian density of a massive vectorial field A^ 

C v ^ l -mlA^-\F^, (7) 

where m v is the mass of the vectorial field and = d^Ay — d v A^ 

(d) the Lagrangian density that corresponds to the Yukawa interaction between the fermionic and the scalar fields 

C Y = -AV^V, (8) 
with A representing the coupling constant of the Yukawa potential. 
Hence, the action of the model in its explicit form reads 



In the above action R denotes the curvature scalar and it was consider 8irG = 1. 

In order to study the evolution of a homogeneous and isotropic spatially flat universe, we use the Robertson- 
Walker metric 

ds 2 = dt 2 - a{t) 2 {dx 2 + dy 2 + dz 2 ) 7 (10) 

where a(t) is the cosmic scaling factor. In this metric, the components of the tetrad, Dirac-Pauli matrices and spin 
connection become 



^ - ^ <='' - — !— , 



s?, r° = 7 °, (11) 
1 



a(t) 



O = 0, ^ = -d(t) 7 l 7°, (12) 
a(t) 2 



where the dots denote time derivatives. 

Here we shall investigate the case where the vector field is time-like, namely, 



A. = (A (t), 0,0,0), (13) 



which is the only possible ansatz consistent with a homogeneous and isotropic universe that leads to a diagonal 
stress-energy tensor with components Tn = T22 = T33. This hypothesis implies that the antisymmetric tensor 
vanishes, i.e., = 0. Furthermore, we shall adopt that the self interaction potential of the fermionic field is given 
by V(ipip) = £ (ipi/i) , where £ and n are constants. The chiral symmetry breaking term reduces to a fermionic 
mass term when n = 2. 

Due to the hypothesis of homogeneity and isotropy the fermionic and scalar fields depend only on time, so that 
we may obtain through a partial integration of ([9]) the point-like Lagrangian: 



C = 3ad 2 — a 3 — 



tpj°ip - tpj°tp + 2iqA tp'-f ip 



a 3 



+ WW ~ \ft + \m 2 b 4> 2 - \m 2 v Al 



(14) 



From the Eulcr-Lagrangc equations for the fields ip and ip applied to the point-like Lagrangian (TT4")) follows the 
Dirac equations for the spinor field and its adjoint, namely, 

ip + \-ip + 1 [qA + nf^r-S + ^ = 0, (15) 
2 a 

^+|-^- ? ^[qA o + <(^)^ 1 7° + A07 ] =0, (16) 
The Klein-Gordon equation is obtained from the Euler-Lagrange equation for <fi and reads 

4> + 3-<p + \^ + mt<j) = 0, (17) 
a 

The component A of the vector field can be determined from the Euler-Lagrange equation, yielding, 

Aq = q 5-, ( 18 

TTly 

which shows that it can be determined once we know the time evolution of ipj°ip. 

By imposing that the energy function associated with the point-like Lagrangian vanishes, namely, 

DC . dC ■ dC ■ „ „ ._„. 

£ = — a + -^0 + ^ip + ip^ - C = 0. (19) 

da d<p dip d$ 



it follows the Friedmann equation 

/ " \ 2 

= p, (20) 



k a > 

where the energy density of the sources of the gravitational field p is given by 



p = ^ + im 2 2 + c(^y i + XipH + \^ (V>7V) 2 - 



Finally, the acceleration equation is obtained from the Euler-Lagrange equation for a: 

a 1 



Ap + Zp). (22) 
a D 



In the above equation the pressure of the sources of the gravitational field reads 



P =-4> 2 - -mi<p 2 + ^ - i)(#r + ■ ( 23 ) 

If we multiply the Dirac equation (fl"5|) by ip and sum with (fl~6|) multiplied by ip we get the following equation 
for the bilinear ipip 



a— 



tpip + tpip + 3-ipip = 0, (24) 



a 



which furnishes through integration 



— C 

tpip = — , where C — constant. (25) 



Hence, the bilinear decays like a prcssureless matter field. 



Instead of using the Friedmann equation we shall deal with the conservation law of the energy density, 
namely, 

p + 3-( P + p)=0, (26) 
a 

which follows from (J20J) and ([22]). 

Now if we substitute the definitions of the energy density density ([2Tj) and of the pressure (|23|) into (|26[) and use 
the Klein- Gordon equation (|17j) we obtain the following differential equation for ^7° tj; 

ijrf ip + + 3- (?7» = 0, (27) 
a v ' 

whose integration leads to 

— C 

= — , where C 1 — constant. (28) 

Hence, we infer from (fT5)) and (J2E]) that the component A of the vector field decays also as a pressureless matter 
field. 

With the results ([25]) and (|28|) we may write the Klein-Gordon (|T7f and acceleration (l22l) equations as 



+ + C ^ = 0, (29) 



a _ 1 

a 6 



a a 

C i<t> & , C 3 , C 4 (3n-2) 



20 2 + ^-C 2 2 + ^ 



(30) 



respectively, where we have introduced new constants 



2q 2 C" 2 

Ci = AC, C2 = m b , C3 = 5 — , C4 = £C 



Furthermore the energy density (|2Tj) and the pressure ([23]) become 



1 • Ci<f> C 2 <\? C 3 C 4 

We note that the pressure equation (|32[) can be decomposed into partial pressures of the scalar p$, fermionic 
and vector pa fields as follows: 

p* = 2^~ _ 2 - ' ^ = ("- 1 )^r> PA = i^- (33) 

The above decomposition will be useful for the identification of the positive and negative contributions related to 
the acceleration field. 

Equations (|29|) and ([30]) represent a system of coupled nonlinear differential equations for the fields a{t) and 
4>(t). A numerical solution of this system of equations can be found by specifying the initial conditions for the 
cosmic scaling factor a(0) and for the scalar field <^>(0) as well as for their derivatives d(0) and 4>(0). 



3 Cosmological results 

We have chosen for the cosmic scale factor a(0) = 1 and for the energy density p(0) = 1. This set of values, combined 
with the Friedmann equation (|2U)) . implies that &(0) = l/v3- Let us assume that at t = the scalar field and its 
time variation are small, so that we choose 0(0) = 10~ 4 and 0(0) = 10~ 2 . These previous choices imply that the 
coupling constant C4 (say) associated with the self interaction potential is determined by equation for the energy 
density (f2"Tj) in terms of the constant C\ , C 2 and C3 , namely, 

C 4 = _ M&f _ d0(O) _ C 2 0(O) 2 _ c 3 
a(0) 3n P[ ' 2 a(0) 3 2 4a(0) (r 

There still remains to specify the free parameters, the exponent n and the constants C\, C 2 and C3. The 
following values have been chosen: n = 1/2, and C 2 = 10 -8 , and we have investigated the role of the vector field 
and of the Yukawa potential in the solutions of ((29)) and ([30]) by varying the constants C\ and C3. 




Figure 4: Scalar field </> versus time t. 



The above parameters are chosen to represent qualitatively the transitions between accelerated and decelerated 
regimes, since we are working with normalized quantities (scale factor, time, energy densities). The system of 
equations that emerge from the original dynamics is highly non-linear, so the main point here is to identify the 
range of values that permit the end of the inflationary period and the beginning of the matter decelerated era, 
taking now into account the presence of the new Walecka-type interaction. 

Numerical integrations show that the cosmic scaling factor leads to an ever expanding universe. The division 
of the cosmological eras can be done in terms of the acceleration field a. Figure 1 shows the solutions for this field 
according to different values of the constants: C± which is proportional to the coupling constant of the Yukawa 
potential and C3, which is proportional to the square of the coupling constant between the fermionic and vector 
fields and inversely proportional to the mass of the vector field. Therefore we have taken the values C\ = 2.0; 2.5 
and C3 = 0.1; 0.2. It is possible to note the existence of three cosmological eras: beginning with an accelerated 
expansion that can be identified with the exit of an inflationary regime, there follows a deceleration and later, a 
dark-energy dominated era when the universe accelerates again. It can be seen that the amplitude of the initial 
acceleration as well as the time length of the decelerated regime are dependent upon the choice of the constants C\ 
and C3: the greater C\ or C3, the smaller the amplitude of the initial acceleration and the longer the deceleration 
of the universe and the shorter initial acceleration. Therefore, the greater the vector field mass and the Yukawa 
coupling the longer the initial acceleration and the shorter the decelerated regime. It is worth noting that in the 
future the deceleration parameter tends asymptotically to a constant value. The search for the values of the free 
parameters and initial conditions that lead to regimes where a transition accelerated-decelerated-accelerated occurs 
is a very hard task, due to the instability of the non- linear coupled system of differential equations (f2T)| and (f5U|) . 

The transition accelerated-decelerated-accelerated can be better understood when we plot the pressures of the 
fermionic, scalar and vector fields (|33p and the total energy density (f3~Tj) as functions of time. These plots are 
shown in Figures 2 and 3. Whereas the amplitude and length of cosmological regimes are sensitive to different 
values of the constant C3, the total energy density and the pressures are not, but only for different values of C\. 
First we note that the pressure of the fermionic field is always negative, while the scalar and vector pressure fields 
are positive. Although the total pressure is always negative, in order to understand the accelerated-decelerated- 
accelerated transition we have to analyze the acceleration field which is given by a = — [p + 3(p0 + P4, + pa)]/6- At 
early times the scalar and vector pressure fields arc small in comparison with the fermionic pressure field so that 
the modulus \p^,\ > (p + Sp^ + 3pa)/3, which leads to a positive acceleration. At intermediate times, the density 
and the modulus of the fermionic pressure decay with time, while the scalar and vector pressures increase. Hence, 
there exists a time interval where \p^\ < (p + 3p^ + 3pyi)/3 and the universe enters into a decelerated regime. After 
this interval the energy density and the pressures of the scalar and vector field become small and a situation where 
\Pip\ > (p + 3j»0 + 3pyi)/3 is recovered, i.e., the universe returns to an accelerated phase. Note that due to the small 
value of C2 the term — C2</> 2 /2 does not contribute significantly to the pressure of the scalar field. The change of 
Ptf, for different values of C\ and C3 is due to behavior of the term </> 2 /2, which can be understood by observing the 
slopes of the graphics 4> X t in Figure 4. Furthermore, there is no significant changes in the pressure of the vector 
field for different values of C\ and C3 chosen here, so that we have represented in Figure 2 only one curve. 

In Figure 4 it is plotted the scalar field as function of time for different choices of C\ and C3 . The decay of <fi 
with time is more accentuated for large values of C\ and C3. In all cases the scalar field tends asymptotically to a 
finite value for large times. The asymptotic behavior of the scalar field (f> for large values of time can be understood 
through the analysis of the Klein-Gordon equation (|29l) . For large values of time the cosmic scaling factor is also 
large so that the Klein-Gordon equation reduces to <fi + 3a(f> / 'a = 0, since C2 was considered as a small quantity and 
the term Ci/a 3 can be neglected. Hence, the integration of the Klein-Gordon equation furnishes cb oc 1/a 3 , which 
is also a small quantity for large times. As a consequence, the pressure of the scalar field becomes very small for 
large values of time and the total pressure is dominated by the pressure of the fermionic field, which is negative 
and contributes to a positive acceleration. It is worthwhile noting that since = identically, the massive vector 
field A v is not free to propagate. 

4 Conclusions 

To sum up, we have proposed a model in which scalar, fermionic, massive vector fields and their interactions 
account for the dynamics and evolution of different cosmological regimes in a homogeneous and isotropic spatially 
flat universe. By observing the behavior of the pressure of the fermionic field with respect to the scalar and vector 
pressure fields, we may say that it is the responsible for the two accelerated regimes, in the early and in the late 
periods of the Universe. As final remarks, we can consider the above investigations as a primer for future work in 
fermionic cosmologies, focusing on the form of the interactions and on testing these spinorial sources using other 
gravity theories. The future possibilities include Walecka-type interactions in scalar-tensor (Brans-Dickc) gravity, 
Bianchi metric fermionic scenarios and supersymmctric inflationary regimes, which is a work in progress. 
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